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ABSTRACT 


Some  old  polynomials  of  L.  J.  Rogers  are  orthogonal.  Their  weight 
function  is  given.  The  connection  coefficient  problem,  which  Rogers  solved  by 
guessing  the  formula  and  proving  it  by  induction,  is  derived  in  a natural 
way  and  some  other  formulas  are  obtained.  These  polynomials  generalize 
zonal  spherical  harmonics  on  spheres  and  include  as  special  cases  polynomials 
that  are  spherical  functions  on  rank  one  spaces  over  reductive  p-adic  groups. 

A limiting  case  contains  some  Jacobi  polynomials  studied  by  Hylleraas  that 
arose  in  work  on  the  Yukawa  potential. 
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SIGNIFICANCE  AND  EXPLANATION 


* 


Spherical  harmonics  are  used  to  solve  many  physical  problems,  especially 
in  potential  theory.  A generalization  of  zonal  spherical  harmonics  was 
introduced  by  L.  J.  Rogers  in  1895.  He  obtained  many  properties  of  these 
polynomials,  including  some  that  would  not  be  found  for  the  classical  spherical 
harmonics  for  another  twenty-five  years.  However  he  was  unaware  that  his 
polynomials  were  orthogonal.  The  orthogonality  relation  is  proved  and 
used  to  obtain  further  results  for  these  polynomials.  Another  limiting  case 
gives  a relatively  recent  result  of  Hylleraas  that  arose  in  his  study  of  the 
Yukawa  potential. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC , and  not  with  the  authors  of  this  report. 


A GENERALIZATION  OF  ULTRASPHERICAL  POLYNOMIALS 


Richard  Askey  and  Mourad  E.-H.  Ismail 
To  the  memory  of  Paul  Turan,  with  respect  and  affection. 


Introduction.  Fejer  [21]  introduced  the  following  class  of  polynomials.  Let 


f(z)  = l 


be  a function  analytic  in  a neighborhood  of  the  origin  with  a real.  Form 

n 

t1-2)  |f(re10)|2  = l P (cos  6>rn  , 

n=0  n 

so  that 

n 

(1-3)  P (cos  0)  = l a a ,cos(n  - 2k)0  . 

" k=0  k n’k 

^n^x)  is  a polynomial  of  degree  n,  and  is  called  a generalized  Legendre  polynomial 
-1/2 

since  f(z)  - (1  - z)  gives  the  Legendre  polynomials.  More  generally, 

f(z)  = (1  - z)  * gives  the  ultraspherical  polynomials  C*(x) . The  orthogonality  rela- 
tion for  the  ultraspherical  polynomials  when  1 > - ~ is 


U.4)  J C*(x)C*(x)  (1  - xV"1/2dx  = 0,  m*n, 

<2X)nr(|)r(X  + i-) 

= n!  (n  + X ) T (X ) ' m = n • 

Fejer  and  Szego  obtained  a number  of  interesting  facts  about  these  generalized  Legendre 
polynomials,  some  of  which  are  summarized  in  [40,  Chapter  VI].  Feldheim  [23]  and 
Lancevickil  [32]  determined  when  the  generalized  Legendre  polynomials  are  orthogonal 
by  showing  that  the  polynomials  must  satisfy  a specific  recurrence  relation.  However 
they  did  not  obtain  an  explicit  representation  for  the  polynomials  and  they  were 
unable  to  find  the  weight  function.  We  will  find  both  the  polynomials  and  the  weight 
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function.  These  polynomials  are  not  new.  They  were  studied  extensively  by  Rogers  in 
the  third  of  an  important  series  of  papers  [36] . The  only  results  from  this  series 
of  papers  that  are  well  known  are  two  identities  that  Ramanujan  rediscovered  and  are 
now  known  as  the  Rogers-Ramanujan  identities.  However  there  are  many  other  interesting 
identities  (some  are  reproved  in  [6]),  and  the  polynomials  in  the  third  paper  are 
probably  more  interesting  and  important  than  any  other  results  in  Rogers'  papers, 
including  the  Rogers-Ramanu jan  identities.  We  will  summarize  some  of  the  results  Rogers 
obtained  for  these  polynomials  and  add  a few  new  ones. 

The  key  to  obtaining  explicit  formulas  is  to  use  the  q-binomial  theorem  and  basic 
hypergeometric  series.  The  binomial  theorem  can  be  written  as 

“ (a) 

(1-5)  (1  - x)_a  = l — ^ xn 

n=0  n! 


(a)  = f(n  + a)/r(a)  = a(a  + 1)  •••  (a  + n - 1)  . 


The  q-binomial  theorem  [5,  Th.  2.1]  is 


“ (a;q)  (ax ; q) 

V n n _ °° 

(q:q>  * ~ ' 

n—u  n °° 


I X I < 1,  |q|  < 1 


(a;g)„  = FT  (1  ~ aq11)  , I q I <1  , 

n=0 


(a; q)  = = (1  - a)  (1  - aq)  •••  (1  - aq 

n . n 

<aq  ;q>„ 


More  general  hypergeometric  and  basic  hypergeometric  series  are  given  by 


fa,  , . . . ,a  i ■»  (a ,)  • • • (a  ) n 

1 r.  _ r 1 n r n x 

r s b , — ,b  ' XJ  (b  ) •••  (b  ) n! 

^ 1 s 1 n=0  In  s n 


(al ar+l 

r+l^r  b,  , . . . ,b 


l (al;q)n  •“  (ar+l,q)n  xn 
" n=0  (bl;q)n  "•  (br'q)n  (q!q)n 


To  motivate  the  study  of  the  polynomials  implicitly  found  by  Feldheim  and 
Lancevickii  we  remark  that  their  generating  function  will  turn  out  to  be 


(6reifl;q)^(Sre  ^jq) 


(1.12) 


a 


, i9  -i6 

(are  ;q)  (are  ;q)  n=0 


y P (cos  e)rn  . 


which  is  a very  natural  analogue  of  the  classical  generating  function 

oo 

<3-33)  (1  - 2r  cos  0 + r2)  * = V C^fcos  0)rn  . 

n 


n=0 


For  1 - 2r  cos  0 + r - (1  - re  ) (1  - re  *9)  , and  the  analogy  between  the  binomial 
theorem  as  given  in  (1.5)  and  the  q-binomial  theorem  in  (1.7)  shows  that  the  function 
(6rel9;q)oo/(arel9;q)  is  a natural  substitute  for  (1  - rel9)’\ 


f 


2.  The  orthogonal  generalized  Legendre  polynomials.  If  a set  of  polynomials  is 

orthogonal  and  satisfies  P ( — x ) = (-l)np  (x) , then  it  must  satisfy  the  three  term 

n n 

recurrence  relation 


(2.1)  2b  xP  (x)  = P (x)  + X P (x),  n = 0,1,...,  P <x)  =1,  P , (x)  = 0 . 
n n n+1  n n-1  0 -1 

Feldheim  [23]  determined  b and  1 for  the  generalized  Legendre  polynomials  that 

n n 

are  orthogonal.  His  result  is 

l _ K . ,K  , > sinh(n  - 1);  _ , 

b - b ^ vb  d ) , . . . . / n - 1 * • • • f 

n 1 10  sinh (n  + 1) £ 

. _ K2  , sinh  (n  - 1)C  . . .2  sinh(n  - 3)f 

n 1 “ 1 1 0 sinh(n  + 1)£  1 0 sinh(n  + 1)£ 

The  coefficients  in  the  original  power  series  for  f (z)  satisfy 


(2.2) 


b = a /a  , 
n n n-1 


Setting  q = exp(-2£),  we  rewrite  these  as 


(2.3) 


(1  n-1. 

b = b,  + <b  - b ) q ■■  ~ S -2 
n 1 1 0 ,,  n+1. 

(1  - q ) 


(2.4) 


An - bi  + 2bi(bi  - bo)q  ;; " 4+i; 

(i  - q ) 


. . .2  2 (i  - q ) 

(bi ' V q 

1 ° (1  - q ) 


With  these  choices  of  b and  X , formula  (2.1)  can  be  written  as 

n n 


2x[b  (1  + q)  - b q - (b  q + b - b )qn]P  (x) 

i u 1 l o n 

= (1  - cin+1)pn+i(x)  * Hb1  * b^q  - bQq) 2 - q"  1 (bj  + b^  - bQ)  2]  P^  (x)  • 

Set  b^  + b^q  - bQq  = a and  b^q  + b^  - bQ  = 8.  This  gives 

(2.5)  2x [a  - 6qn]P  (x)  = (1  - q"+1)P  (x)  + [a2  - 02qn_1]P  , (x)  . 

n n+i  n-i 

To  find  the  polynomials  P^fx),  multiply  (2.5)  by  rn+1  and  sum,  recalling  that 

P_^ (x)  = 0,  and  PQ(x)  = 1.  If 

oo 

f(r,x)  = l P (x)r"  , 
n=0  n 


i 


the  resulting  equation  is 


2xarf(r,x)  - 2x6rf(qr.x)  = f(r,x)  - 1 - (ffqr.x)  - 1]  + a2r2f(r,x)  - 62r2f(gr,x) 


(2.6) 

Set  x = cos  0 and  rewrite  (2.6)  as 


2 2 2 7 

(1  - 2xar  + a r )f(r,x)  = (1  - 2xBr  + 6 r )f(qr,x) 


(2.7) 


f (r  ,cos  6)  = i-1--8™  )(1  - 6re_^>  f (gr ,cos  6)  . 

(1  - are  ) (1  - are  ”) 


Iterate  (2.7)  to  obtain 


(2.8) 


(Brel9;q)  (fire  i6;q) 

f (r.cos  0)  = 2 _ £1  f(qnr(COS  0)  _ 


(arei0;q)  (are  i0;q) 
n n 


Then  let  n -*•  The  result  is 


(2.9) 


f (r ,cos  0)  = 


(Bre10;q)  (Bre  i6;q)  •» 


. i9  . . -i0  u n 

(are  ;q)  (are  ;q)  n=0 


= \ P (cos  0)rn  . 


This  is  the  generating  function  stated  in  (1.12). 

The  above  derivation  was  formal,  but  it  is  easy  to  justify  it,  since  the  function 
f(z,cos  0)  is  analytic  for  |z|  < |a|  A.  There  is  no  loss  in  generality  in  taking 
a = ^or  it  can  be  removed  by  scaling,  with  this  simplification  define  C (x;6|q)  by 


(2.10) 


(Bei9r;q)oo(Be  iSr:q)o 
n=0  (ei9r;q)  (e  i8r;q) 


I C (x;6|q)rn 


We  will  use  x and  cos  0 interchangeably.  When  6 = qA  it  is  easy  to  see  that 


(2.11) 


lim  C <x;qA|q)  = C* (x) 
q»l"  " n 


This  follows  from  either  the  generating  function  (2.10) , using 
\ . \ 


(2,12) 


lim 

<T*l”  (x;  q) 


(q  x; q)  « (q";q)  „ (*) 

= = xim_  \ n _ t n ~n 


#4 / , - " 

00  q-*l  n=0 


(q;q> 


I 

n=0 


nl 


(1  - x) 


or  from  the  recurrence  relation 

(2.13)  2x(l  - Bqn]Cn(x;8|q)  = (1  - qn+1) Cn+1 (x ; B | q)  + (1  - bV*'1)^^  (x;  B |q) 


-5- 


which  becomes 


XX  X 

(2.14)  2x(n  + X)C  (x)  = (n  + 1)C  (x)  + (n  + 2X  - 1)C  _ (x) 

n n+1  n-1 

when  6 = q\  every  term  is  divided  by  1 - q and  q -*•  1 . 

If  C (x ; 6 1 q)  = k r (x;8|q)  with  k = (8;q)  /(q,q)  , then  (2.13)  becomes 
n n n n n n 

, , ft2n— 1...  n. 

(2.15)  2xr  (x)  = r (x)  + - — ~ - 3 1 ~ q r (x)  . 

n n+1  (1  - Bq"-1) (1  - 6qn)  n_1 

Now  we  see  that  the  case  |q|  >1  can  be  reduced  to  the  case  |q|  <1  when  8 is 
replaced  by  8 *.  Thus  there  was  no  loss  in  generality  in  assuming  |q|  < 1.  We  use  6 

rather  than  q^  since  8 could  be  negative,  and  the  case  8=0  is  very  important. 

The  conditions  for  orthogonality  can  be  obtained  from  (2.15).  See  [20],  [27,  Ch.  II,  Th.  1.5]. 
It  is  necessary  and  sufficient  that 

a - fV;1!. ix. - > o,  

(1  - 6q  > (1  - 8q  ) 

-1/2 

When  0 < q < 1,  the  cases  n = 1 and  n = 2 imply  -1  < 8 < q , and  these 
conditions  are  easily  seen  to  be  sufficient.  When  -1  < q < 0 the  cases  n = 1,  n = 2, 

and  n = 3 imply  -1  < 6 < -q  ^ and  these  conditions  are  sufficient.  The  trivial 


case  q = 0 holds  when  8 > -1.  The  only  way  this  inequality  can  hold  when  q -*•  1 
is  to  have  8 = q* , in  which  case  the  condition  is  X > - ^ . Finally  the  case  q -<•  -1 
can  only  hold  when  8 = |q|^  and  X > -1  or  6 = ~|q|  and  X > 0. 


Explicit  representations  for  the  continuous  q-ultraspher ical  polynomials.  There 
is  another  set  of  orthogonal  polynomials  that  are  basic  hypergeometric  series  and  have 
the  ultraspherical  polynomials  as  limits  when  q 1.  These  will  be  called  the 
discrete  q-ultraspher ical  polynomials,  because  their  distribution  function  is  a discrete 
measure.  See  [7],  The  weight  functions  for  the  polynomials  under  consideration  in 
this  paper  are  absolutely  continuous  when  1 3 1 <1,  so  these  polynomials  will  be 
called  the  continuous  q-ultraspher ical  polynomials.  The  adjective  will  be  used  in 
both  cases,  since  it  is  not  clear  which  will  be  more  important. 

To  find  an  explicit  representation  use  the  q-binomial  theorem  (1.7)  in  the  generat- 
ing function  (2.10).  The  result  is  (x  = cos  6) 


(3.1) 


C (x;6|q)  = l 


n (6;q)k(B;g)n.k  i(n_2k)9 

k=0  (q;q)k,q;q)n-k  6 


n (6;q).  (6;q) 

= 1 r — j — cos(n  - 2k  )0 

k=0  (q;q)k(q;q,n-k 

lB;q)n  me  ? «fn;q>k<e.-q)k  _2ie.k 

= e i i-n  -i <qe  6 > 

4 4 n k=0  (qx  8 ;q>k(q;q)k 


since  (a;q) 


(3.2) 


(a ; q)  q 
n 


k(k+l)/2 


n-k  n.k,  1-n  -1  , 

(-aq  ) (q  a ;q) 


. This  gives 


k 

(6;  q) 

C (cos  0 ; 8 | q)  = 

n ' (q;  q) 


n in9 


. f 

2 1 


~n  , 

q 


l-n-l 
q 8 


„-i  -2ie 
q,q8  e 


Another  representation  can  be  obtained  as  a special  case  of  a more  general  set  of 
orthogonal  polynomials  considered  in  (12).  The  polynomials  are 


(3.3) 


P (x)  = * 

n 4 3 


-n  n-1  i0  -i0 

q , q abed,  ae  , ae 


; q.q 


ab,  ac,  ad 


x = cos  0.  Their  recurrence  relation  is 


(3.4) 


2xP  (x)  = a P (x)  - (a  + c - a - a P (x)  + c P , (x)  , 
n n n+l  n n n n n-1 


-7- 


where 


(3.5) 


(1  - abq")  (1  - acqn)  (1  - adq°)  (1  - abcdg”  *) 
a(l  - abcdq2n  S (1  - abcdq2n) 


(3.6) 


a(l  - qn) (1  - bcgn  1 ) ( 1 - bdqn  1 ) ( 1 - cdq”  1 ) 
(1  - abcdq2”  2) (1  - abcdq2n  *) 


First  set  c = -a,  d = -b.  In  this  case  the  weight  function  is  even  and  a + c = a + a 

n n 

1/2 

Then  set  b = aq  . With  these  specializations,  the  coefficients  a^  and  c^  are 


-] 


(3.7) 


,,  4 n, 

3 = ±L-~  a <1  > 

n a(l  - a2qn) 


(3.8) 


= a(l  - qn) 

, , 2 n, 

(1  - a q ) 


n 2 4 

Set  P (x)  = k s (x)  with  k = a (a  ;q)  /(a  ;q)  . The  recurrence  relation  (3.4)  becomes 
n n n n n n 


(3.9) 


2xs  (x)  = s (x) 
n n+1 


4 n-1,  n, 

4-  **  , ’ (1  ~ * > S (X) 

(1  - a2q"  X) (1  - a2qn)  n 1 


This  gives  r (x)  = s (x)  with  8 = a , where  r (x)  satisfies  (2.15).  Combining 
n n n 

the  above  formulas  we  see  that 


(S2;q) 

(3.10)  C (cos  0;  6 |q)  = c 51 

n 8 ' (q;q> 


4*3 


q'n,  qV.  01/2ei6,  81/2e-i6 

6q1/2,  -8g1/2,  -8 


; q.q 


Another  way  of  writing  (3.10)  is 


k-1 


(3.11)  C (x;  8 1 q) 


(62;q) 


<q  n;  q)  v (qn82;q)  ] f (1  - 261^2xq;’  + 6q 


23, 


l 


k 

j.-P. 


8 ' “ (q;  q)  k=0 
n 


(Bq1/2.q)k(-Bq1/2;q)k(-8;q)1.(q;q)k 


This  is  probably  as  close  as  one  can  get  to  a single  sum  that  gives  these  polynomials 

as  a series  in  the  polynomial  variable  x Another  formula  which  shows  the 

polynomial  character  of  C (x;B|q)  is 

n 


(3.12) 


n (Bsq).  (8;q)  . 

C (x;B|q>  * J 7 r~7 T 

n k®0  <9»«Dk<q,q,n-k  n_2k 


(x) 


where  Tn(x)  is  the  Chebychev  polynomial  of  the  first  kind  defined  by  Telcos  8)  = cos  n6 


-8- 


Two  interesting  formulas  follow  from  (3.12).  One  is  the  special  case  6 = q. 


(3.13) 


The  second  is 


C (cos  9 ; q | q ) = £ costn  - 2k)  9 = sin(n  V>9  ■ 

n sin  9 


(3.14) 


C (cos  9 ; 8 ! q) 

, . n.  n 

lim  (1  - q ) — — — = COS  n9 , 


B+l 


2(1  - 6) 


n = 1,2,...  . 


Observe  that  in  both  cases  the  polynomials  are  independent  of  q. 

In  the  case  of  ultraspherical  polynomials  the  zeros  of  the  polynomials  C*(x) 

n 

have  absolute  values  that  decrease  as  X increases.  This  suggests  that  the  zeros 
of  c^tcos  9;q*|q)  lie  between  the  zeros  of  cos  n9  and  sin(n  +1)9  when  0 < X < 1 
and  0 < q < 1.  Feldheim  (23)  has  shown  this,  but  his  notation  is  so  much  at  variance 
with  ours  that  we  give  a new  proof.  By  a theorem  of  Fejer  (22),  [40,  Theorem  6.5.2] 
it  is  sufficient  to  show  that  the  coefficients  are  a moment  sequence  on  !b,l],  i.e. 

1 

a = / tndu(t),  dp(t)  > 0 . 


Thus  we  need 


(3.15) 


(q*;q)  1 


Let  dp  have  a point  mass  equal  to 


.X  1-X  n+X , 

(q  .-q),_<q  ;q>„  t <tq,-q) 


...  1_x  > 
(tq  ; q) 


at  t = q , 


0,1,...  . 


When  0 < X < 1 this  mass  is  positive,  and  the  q-binomial  theorem  (1.7)  shows  that 
(3.15)  holds.  The  inequalities  for  the  zeros  are 


*(k ' V X 

£_  < eA 


kn 


k,n  = n + 1 ' 


(3.16) 


where  9 


0 < X < 1,  0 < q < 1,  k - 1 JjfJ  , 


< 6X 

n + 1 = k,n 


"(k  - j) 


0 < X < 1,  0 < q < 1,  1 


III- 


X 

k,n 


are  the  zeros  of  C (cos  6 : I q)  ordered  by  0,*  < 0,*  . There  is 

n ■ k,n  k+l,n 


equality  in  (3.16)  only  when  n is  odd  and  k = (n  + l)/2. 


-S- 


An  unlikely  looking  formula  follows  on  equating  the  basic  hyper geometric  series 


in  (3.2)  and  (3.10).  It  can  be  written  as 


-n  2 
q , a 


„ ; q,  S5_  = 

2 1 1-n  -2  2 

q a a 


2 \ xn(a4;q) 


r -n  n 4 ”1 

I q , q a , ax,  ax 


n . 2 x 4V3  2 1/2  2 1/2  2' 

a (a  ;q)  (a  q , -a  q , 'a 


When  a = q and  the  limit  q •+  1 is  taken  the  resulting  identity  is 


-n,  a 

xn  (2  cx) 

-n,  n + 2a  , , .2 

2 

n F 

(1  - x) 

; x 

(a)  2 1 

1 ' 4x 

1 - n - a 

n 

“ + 2 

This  is  one  of  the  iterated  quadratic  transformations.  The  first  quadratic  type 
transformation  for  basic  hypergeometric  series  seems  to  be  that  of  Carlitz  [17] . 
Another  is  given  in  (7],  where  the  discrete  q-ultraspher ical  polynomials  are  related 
to  some  discrete  q-Jacobi  polynomials.  Another  will  be  given  in  [12]. 

An  interesting  inequality  follows  from  (3.1), 

(3.19)  |c  (x;B|q) | < c (l;B|q) . -1  < B < 1,  -i  < q < 1 • 

For  | cos  (n  - 2k.)  0 1 < 1 and  (B;q)k  > 0,  (q;q)k  > 0 when  -1  < B < 1/  ”1  < q < 1- 

Unlike  the  classical  case  of  C*  (x) , when 

n 


it  is  impossible  to  find  the  value  C (l;B|q)  as  a simple  product.  There  are  two 
interesting  points  where  the  value  can  be  given  as  a product.  From  (2.10) 

. n (Bir ; q) ^(-Bir;  q)  ^ (-S2r2,g2)(c  - (fiW)  n _ % n 2n 

l (0;  6 1 q)  r - (ir;q)  (_ir;q)  2 2.  ^-.2  2.  1 


(-r  ;q  > n= 0 (q  ;q  ) 


(S2.-q2) 

C (0;S|q)  = (-1)"  ■- 

20  (q'  :q  )n 


c2n+i(0;8|q)  = 0 . 

a,  b 

Heine  [28],  see  also  [5,  Cor.  2.4]  found  an  analogue  of  Gauss'  sum  of  s 1 


I c 


l^;<J  (C^>o 


If  the  series  terminates  this  is  the  correct  value  without  the  condition  - 

ab 


< 1. 


Set  e 


2i6 


=6  in  (3.2).  This  gives 


(3.23) 


81/2  ♦ B-1/2 


-•  P I q 


(62,-q) 


(q;q) . 


n B-n/2, 


0 < 8 < q 


-1/2 


When  1 < 0 < q 


1/2 


this  is  interesting,  for  then 


lim  C 

r 

n-+«) 


3 1/2  + B~1/2 


; B|q 


= 0 . 


This  goes  to  zero  fast  enough  so  the  distribution  function  has  a point  mass  at 


(8  +6  )/2  when  1 < 3 < q ^ ~ , see  (33)  . By  symmetry  there  is  also  one  at 


x = “(6  ^ + B " ) /2 . These  are  the  only  point  masses  when  1 < 3 < q 3^2. 


For  convenience  the  first  four  of  these  polynomials  are  given  next. 
C0(x;B|q)  = 1 

2(1  - 3) 


Cj (x ; 6 1 q)  = 


(1  - q) 


(3.24)  C2(x;8|q) 


= 4(1  - 3)  (1  - 8q)  2 _ (1  - 


“) 


(1  - q)  (1  - q ) (1  - q ) 

2, 


C (x;  8 1 q)  = 8(1  - 3)  (1  - 8q)  (1  - Bq‘)  x3  _ 2x  ( 1 - 3)  1 2 + q + B (1  - q*)  -Bq(l  + 2q)]  _ 
(1  - q)  (1  - q2)  (1  - q3)  (l-q2)(l-q3) 


From  this  point  on  the  formulas  become  more  complicated. 


Another  useful  expression  for  (cos  0;8|q)  can  be  obtained  by  applying  the 


q-analogue  of  the  Pfaff-Kummer  transformation  to  (3.2).  This  transformation  is 


(3.25) 


't 

2 1 


a,  b 


: q.x 


( ax ; q)  - [b  ' q)  (a; q> n (-*b) " 

* ao  ^ V J n 


<2> 


(x;  q) 


«>  n=0 


(c;q)  (q;q)  (ax;q) 
n n n 


See  Andrews  (4)  for  a proof  of  this  formula.  The  resulting  formula  is 


(3.26) 


(6;  q) 

C (cos  0 ; 3 1 q)  = r 

n (q;q) 


n inO  -1  -2i6  1-n  . 

- e (3  e q ; q)  * 


(^)  k 

, -n  . , 1-n. -2  . '2'  , -2i0 

n (q  ;q)k(q  6 ;q)kq  (-qe  ) 

. „ , l-n.-l  , , 1-n  -1  -2i0 


k=o  (q  3 ;q)k<q  8 e ; q) k (q» q) k 

Now  sum  this  series  in  the  opposite  direction,  i.e.  replace  k by  n - k,  and  simplify 
to  obtain 

2i0 


(3.27) 


,„2  . -ine 

(3  ;q)  e 

Cn(cos  0 ; 6 | q)  - n-  ■ ~ - 3*>2 


-n  . . - 

q . B.  Be 


(q;q)nB 

From  this  it  is  easy  to  obtain  the  generating  function 


B2,  0 


> q.  q 


” C (cos  0 ; 6 | q)  (,) 

(3.28)  £ q (8r)  = (-re  ;q)m 

n=0  (6  ;q) 


'»  2 1 


B.  Be 


2i0 


; q,  -re 


-i0 


1 


4.  The  orthogonality  relation.  The  orthogonality  relation  for  the  q-Wilson  polynomials 
defined  in  (3.3)  is 

00 

P (x)P  (x)  TT  (1  - 2(2x2  - 1) qk  + o2k) 
l n m , _ 

(4.1)  / — — - dx 

. h(x,a)h(x,b)h(x,c)h(x,d)  / 2 


2ira2ribc;q)n(bd;q)n(cd;q)n(q;q)n(l  - abcdq"”1) 

= • (ab;q)n(ac;q)n(ad;q)n(abcd;q)n(l  - abcdq2"'1) 

(abcd;q) 

00 

(ab;q)  (ac;q)  (ad;q)  (bc;q)  (bd;q)  (cd;q)  (q;q) 

<*>  00  00  OO  00  oo  c 

when  |a|,|b|,|c|,|d|  < 1,  where 


h(x,a)  = "]  f"  (1  - 2axqn  + a2q2n) 


In  general  there  will  be  finitely  many  mass  points  outside  (-1,1).  See  [121  for  details. 

For  the  continuous  q-ultraspherical  polynomials  when  |s|  < 1,  the  specializa- 
tions of  the  last  section  lead  to 


K 


(4.3)  / C (x;6|q)C  (x;8|q)w  (x)(l-xV1/2dx  = 2t  |-B  ' ^ ^ “ s 

-1  m 6 (1  - 6q")  ‘^’n  (62;q)  (q;q,  " 


w<x)  » fr  -i1-.  2<2*2  - * s2k> 

8 7 k 7V  7 

k=0  (1  - 2 (2xZ  - l)q*B  + q 6 ) 


This  can  also  be  written  as 


. 2i6  . . -2i0  , 

(e  ;q)Je  ,-q)^ 

w.(cos  0)  * — — — — 

B 2i0  -2i0 

(Be  ;q)  (Be  ;q) 


To  get  a better  idea  of  this  weight  function  observe  that 

u -1/2  _ ,2,  . .2 


W (x)  (1  - x ) = (4  - 4x  q ( (1  + q)  - 4qx  ) (1  - xZ) 


When  q » 1 this  becomes  42(1  - x)2)2'1/2. 


2. -1/2 
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I 


Before  considering  other  special  cases  we  will  give  a direct  proof  of  the 
orthogonality  relation  (4.3).  There  are  two  reasons  for  this.  One  is  that  no  other 
proof  has  appeared  yet,  and  the  three  other  proofs  that  we  know  use  results  that  are 
not  needed  in  the  present  proof.  Also  this  proof  will  give  a direct  proof  of  an 
important  identity  of  Rogers  [36). 

As  a first  step  we  compute  a trigonometric  moment  that  is  of  independent  interest, 
/ 2ik6  rv  2ik0  -ry  (1  - r2iPgn)  (1  - e 2l°qn)  .. 

/ e w (cos  0)d0  -J  e IT  ~ 2i0fl  n,  “ -2i0~  n.  d®  * 

0 0 n=«0  (1  - e gq  ) (1  - e 8q  ) 


Use  the  q-binomial  theorem, 


, 2i0  1 /a-1  \ 

(e  ; q)  °»  (6  ;q) 


(Be2i0;q)  n=0  (q;q)n 


n .n  2in0 
— Be 


to  get 


2ik0  . ....  ; <8'1;q)n  .»  r l6'1;q)m  rm  / 2i(k+n-m)0.n 

/ e w (cos  6)d0  = l - B l ■, q , - B J e d0 

0 n=0  4 ^ n m=0  4 4'm  0 


r (S'l!q)n(6'1;q)ktn  .k+2n 
” n=0  (q;q)n(q;q,k+n 


„k  (6"1;q’^  * 

= "6  2*1 


S'1,  qV1 


) q,B 


k+i 


Heine  (281,  [29,  p.  1061,  see  also  (5,  Cor.  2.2],  gave  an  important  transformation  of 
the  general  , , 


a,  b 

(ax;q)  ^ (b;  q)  ^ 

re 

b* x 

(4.6) 

2 1 

; q,x 

c 

(x;q)  ^(c ; q)  ^ 2*1 

; q,b 

ax 

The  iterate  of  this  is 


. . c 

abx  . 

a,  b 

(bx; q) ^ (^  ; q). 

, b 

c c 

(4.7) 

4 

t 1 

? q,x 

c J 

q)  (c-q)„ 

v? 

2 1 

; q'  b 
bx 

Use  (4.7)  on  the  * above  to  get 
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; ,f  '-1  J l* xuPUgpi  I*  ^ -,*-*y  wwjfMWt 


; q,6q 


k+l 


Summing  this  ^ (there  are  only  two  nonzero  terms)  gives 


(4.8) 


2ik0  ”0  (B  1;q)  (X  + qk)  (B;q>  (Bq;q) 

/ e"1"  w Icos  0)d0 ^ = 

n 6 (6q;q)  -2 


(q;q)  (6  ,q) 


|s|  < 1.  The  condition  |8|  <1  was  used  to  expand  w. (cos  0)  by  the  q-binomial 

P 

theorem,  and  the  series  that  was  used  does  not  converge  when  02  > 1.  The  case  6=1 


is  trivial,  since  w^ (cos  0)  = 1.  In  this  case 


I e 
0 


2ik0 


w^ (cos  0)d0  =0,  k = ±1,42,...  , 


n , k =•  0 . 

To  continue  the  proof  of  the  orthogonality  relation  consider 


/ C (cos  0;8|q)T  (cos  0)w  (cos  0)d0 
q n n”*K  p 


where 


T (cos  0)  = cos  n0 
n 


This  integral  is 


n (B;q)  . ( 6;  q)  . n . , 

V 3 n~3  f 1 (n-2j ) 0 

j=0  (<J!q)j(qiq>n-j  0 


i (n-2k) 0 ^ -i (n-2k)0 
e + e 


w (cos  0 ) d0 

p 


= - l (Blq)j-B;q)n-3  f (e2i(n-j-k)0  2i  (k-j)6 

2 jto  0 6 


0)d0 


(6;q)  (8q;q) 


2 (62;q)„(q;q), 


[i  . + i. ) 

n-k  k 


where 


n q)  ^ (P» q)  n_ j B^^B’Sq)^ 


j=0  (qiq)  jtq.-qJn.j  (Bq.-q)^ 


(1  + qk"j) 
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Use 


(a;qjizi  = (a;q) n <q1~">»'1»q)i  fb]j 


(b;q)  . (biq)  1-n  -1  ia 

n-j  n (q  a ,-q)..  1 


to  rewrite  this  as 


<B;q)  (B_1iqj  n '-~n 

(4-9)  I,.  - Bk(H-qit)  l - 


k (q;q)n(6q;q)k 


-k.-l 

J.'4-6  8 -3 


l-k 


k n (q  ; q)  ^ (q  B j q)  (8;q)  ..  (-q  ,-q).. 


j=0  (q.-qj^fq1  °B  1;q)^(q1  kB;q)j(-q  k;qK 


^-q3 


This  sum  is  a balanced  if ^ (i.e.  the  product  of  the  numerator  parameters  times  q is 


the  product  of  the  denominator  parameters  that  are  listed  in  the  and  so  may  be 


transformed  by 


(4.10) 


4*3 


-n  . 

q , a,  b,  c 


i q.q 


d,  e,  f 


de  fdf 

be  ' q be  ' q 

'n'  n 

-n  d d 

q - a-  r < r 

c b 

(e;q)  (f ;q) 
n n 

4*3 

; q.q 

de  df 
' be  ' be 

where  q1  °abc  = def . See  [7]  or  [12]  for  this  transformation.  Let  d = -q  k,  a = B. 


This  gives 


(4.11) 


ia  /o-1  , 0k-n  . k,  (l-k)n,  k-n  , ,.2 

(8;q)  (B  ;q).  B (1  + q )q  (q  ;q)  (B  ;q) 

y a 0 x n n 

k 1-n  -1  , , l-k. 


(q;q)  (Bq;q)  (q  6 ;q)  (q  8;q) 

n k n • 


* 

4 3 


q ",  8,  q 1,  -6 


-k  k-n  2 1 q,q 

-q.q  .6 


(Bq) 


,,  n-2k,  -1  , ..2  , . l-k  , 

k (i  - q ) (B  ;q)k<8  ,-q)n<q  sq> 


n-1 


(1  - q1  kB)(Bq;q)  (q;q)  (q2  kB»q>  . 

K n n-1 


Thus  I = 0,  k = 1,2, ...,n  - 1,  and  I - I = (8  ;q)  /(6q;q)  . Since 


T (x)  = 2n_1xn  + •••,  n - 1,2,...,  T (x)  = 1 , 
n 0 


and 


(B;q) 


Cn -2"<55T»“*  •••■  


it  is  easy  to  check  that  (4.3)  holds. 
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M m 


The  argument  above  is  easy  when  |8|  < 1,  it  would  be  more  complicated  for  the 
values  of  6 with  | 8 | >1  when  the  polynomials  are  orthogonal  with  respect  to  a 
positive  measure.  Rather  than  try  to  carry  out  this  argument  (we  have  not  tried) or 
introduce  the  methods  that  we  know  will  lead  to  this  orthogonality  we  will  now  use 
the  orthogonality  to  find  some  important  identities.  The  complete  orthogonality  rela- 
tions will  be  given  in  til)  and  [12].  Different  methods  are  used  in  these  two  papers. 

The  argument  above  did  not  need  k = 0,1,..., n.  Let  k be  successively  -k 
and  n +■  k in  (4.11).  A calculation  gives 


(4.12)  / C (cos  0 ; 6 I q) T (cos  6)w  (cos  0)d8 

q n n+^K  p 

^6k(8"1.-q)k(q;q)n+k  (1  _ qn+2k}  (8;  q)  „ (6qn+k+1;q)0 

(wqMqjq)  ,,  n+k,  . . ..2  n , 

k n (1  - q ) (q.-qJ^B  q ; q)  ^ 

Using  the  orthogonality  relation  (4.3)  and  (4.12)  gives 


k = 1,2,... 


j . (6-1;q).  (q; q)  . (1  - qn)  (1  - Sqn'2k) 

Vx)  " l 6 — — : - . c„-2k(*'-eh>  • 


k=0  (q;q).  (6q;q)  (1  - q ) 2(1  - 6) 


Then  if  (4.13)  is  used  on  the  right  hand  side  of  (3.11),  the  result  is 
^ (Y;q)k(Y;q)n-k(q,q)n-2k  n-2k 

<4-14)  C"(X”k)  = Jo  ^>qU  U ' * > • 


f_n-2k  n/2-k+l  n/2-k+l  n-k  „-l  -k 

q . q . -q  . rq  , 8 , q 


6*c  * q-  “ 

n/2-k  n/2-k  -1  1-k  n+l-2k  n+l-k  Y 

|q  , -q  , y q , 8q  , q 


8a  Cn-2k(x;Blq) 


(1  - 8) 


The  is  very  well  poised  and  so  can  be  summed  by  a theorem  of  Jackson  [37,  p.  96) . 


The  resulting  series  is 


LfJ 


C (x;  8 1 q)  = l 6 


k(Y8  .q)k(Y;q)n_k 

(q;q)k(Bq;q)n_k 


(1  - 6)  n-2k 


- c ,vtx'-e|q)  • 


This  very  important  formula  was  found  by  Rogers  [36].  He  obtained  it  by  finding 
the  coefficients  for  small  values  of  n,  guessing  the  answer,  and  then  proving  it  by 
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induction.  The  special  case  when  y = 0 and  8=1  was  used  in  his  second  paper  [35] 
to  obtain  the  identities 


n <q.-q  > (q  .-q  > 


“ n +n 

V -3 1 

*•_  <q;q)  ,2  5.  , 3 5, 

n=0  n (q  ;q  ) (q  ;q  ) 


While  it  is  not  particularly  easy  to  read  and  understand  these  papers  of  Rogers, 
we  have  found  it  easier  to  read  them  after  understanding  the  polynomials  of  Rogers. 

It  is  also  easier  to  read  them  in  the  order  [36],  [35],  [34]  rather  than  the  order 
in  which  they  were  written.  These  papers  are  very  interesting  and  contain  many  other 
important  results.  Probably  the  most  interesting  is 

min(m,n) 

(4.18)  C (x;8|q)C  (x;  6 1 q)  = l a(k,m,n)C  (x;  6 1 q) 

m n . _ m+n-zK 


a(k,m,n)  = 


^;^mtn-2k(B;^n-k(S;^m-k(S--q>k(B2^>mtn-k(1  ~ 

(62;q)m+n-2k(qjq)n-k(<3;q)m-k(q;q,k(e<I;q,m+n-ka  * 61 


Rogers  found  this  result  in  the  same  way  he  found  (4.15),  by  working  out  the  coefficients 
for  small  m,  guessing  the  answer,  and  proving  it  by  induction.  We  do  not  have  a 
better  proof  of  (4.18)  at  this  time. 

Rogers  [36]  pointed  out  that  the  special  case  q = 1,  8 = q*  is  an  important 
result  for  spherical  harmonics.  However  no  one  picked  up  this  result  from  [36]  and 
the  special  case  8 = q*,  q = 1 of  (4.18)  was  next  stated  by  Dougall  [19]  almost 
twenty-five  years  later.  The  only  special  cases  that  were  known  before  Rogers  found 
the  general  result  were  the  trivial  cases  6=1,  which  is  equivalent  to 
cos  n0  cos  m0  = ^ [cos(n  + m)0  + cos(n  - m)0]  , 


6 = q,  which  is 


sin(n  + 1)6  sin(m  ♦ 1)8  _ nUn(ym'n>  sin(m  + „ + 1 - 2k)  6 

sin  e sin  6 * sin  0 

k*0 


and  one  nontrivial  case,  B - q ' , q * 1, 


min(m,  n)  m+n  + ^-2k  2,2  ,2  <">  + n - k)  I 

p <x)p  (x)  = i 1 Iv -j»-r _!"-!< 

n m , *•_  1 fiT 

k-0  m + n + - - k kl  (m  - k) I <n  - k) ! | 

^ J m+n-k 


Pm+n-2k (x)  ' 


a result  of  Ferrers  [24,  Example  10,  p.  156J  and  Adams  [1).  See  Chapter  5 of  [9]  for 
a summary  of  some  of  the  known  results  on  the  linearization  problem. 

Formulas  (4.15)  and  (4.18)  can  be  inverted.  The  inverse  of  (4.18)  is 

(8;q)m<B;q)n  *®;<I*m+  B,in(B,'n) 

<4-19)  (q;q)m(q;q)  Cm+n(X!®^)  " TqTq)  ~ l b (k'm'n)  Cm-k  (X>  6 I q)  Cn-k  (X;  6 lq)  , 

m n m+n  k*0 


b(k,m,n) 


. -m-n  -2  , . 2k-m-n„-2,  -1  , _\k 

(q  6 »q> k (1  “ q B ) (6  ;q)k  g2 
(qjq).!!  - q'""V2)(q1_*"V1»q). 


The  limiting  case  8 - q , q - 1 of  (4.19)  was  found  by  Al-Salam  [3] . A proof 

of  (4.19)  can  be  given  following  his  proof;  use  (4.18)  to  replace  C . (x;6|q)C  , (x;8|q) 

m-k  n-k 

and  invert  the  order  of  summation.  The  resulting  series  is  a very  well  poised  , 

6 5 

which  is  summed  by  Jackson's  result.  Details  are  left  to  the  reader. 


The  other  inverse  is 


'6(x,Cn(x;  6^q)  " ^ a(k,n)C^,k(x;y|q)w^(x)  , 


a(k,n) 


Bk(Y/B;q)k(qn+1;q) 2k <Y2q"+2k; q)^ (8qn+k+1; q)^ (6; q) m (1  - Yqn+2k) 
(q»q)k(Yqn+k;q)00(82qn;q)0>(T!q)00 


This  follows  from  (4.15)  by  the  general  argument  given  by  Askey  in  [8] . The  series 
(4.20)  converges  when  |b|  < 1,  |q|  < 1.  For  |a(k,n) | - 0(Bk)  and 


(y*  q)  v 

lCn*2k<X’^q,l  i(nUI  “x  7^ 


A (Y  ,q)  (n  + 1) 


for  Y.q  fixed,  |q|  < 1. 


Rogers  [36J  found  a number  of  other  formulas.  The  reader  is  referred  to  [36]  for 
these  results. 
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5.  Special  cases.  It  is  not  surprising  that  interesting  results  are  found  when  q •+  1 . 
However  it  is  surprising  that  an  interesting  result  could  be  found  by  letting  q ■*  0. 

To  see  this  consider  the  generating  function  (2.10)  when  q = 0.  It  is 

(5.1)  1 - 2Br  cos  6 = (1  Bre1^  (1  - Bre^6)  . J c (cos  9;6.0)rn  . 

1 - 2r  cos  0 + r (1  - re1  ) (1  - re  1 ) n=0  n 


2 sin(n  + 1)9 


n = 1, 2,  . . . , 


C (x;  6 1 0)  = 26(1  - 6)T  (x)  + (1  - 6)‘u  (x)  , n = 1,2,...  , 

n 1 n n 

CQ  (x;  6 1 0)  « 1 . 

These  can  also  be  written  as 

C (cos  0 ; 6 1 0)  = <1  - 6)  Sln(n  + ^-9  - B (1  - B)  Sln(n  ~ 1)9  , n = 1,2 

n sin  v sin  o 

CQ (cos  0; B|0)  =1 


A partial  fraction  decomposition  gives 


C (cos  9 ; 3 0)  = 26(1  - 6) cos  n0  + (1  - 6) 
n 

CQ (cos  0 ; 6 | 0)  = 1 


C (x;  8 0)  = (1  - B)V  lx)  - 8(1  - 6>U  .lx),  n = 1,2,3,...  , 

n n n-2 

CQ (x; 6 | 0)  = 1 . 


The  orthogonality  relation  is 


(5.3)  / C (x;  8 1 0)C  (x;  6 1 0) 


n 2.1/ 2 

— — 1 r dx  = - (1  - 6)  6 , (m,n)  * (0,0) 

(1  t 6) 2 - 4Bx2  2 


1 , 2,1/2 

I (C  (x;6|0))2  — ^ l- Y1 dx 

-1  (1  + 6) ^ - 4Bx^ 


2(1  + 6) 


when  6 < 1 . 


This  is  a known  result.  It  follows  easily  from  the  Szego-Bernstein  theory  [14], 
[40,  52.6];  it  was  also  given  by  Karlin  and  McGregor  [31).  When  6 > 1 there  are 


two  point  masses  that  need  to  be  added  outside  [-1,1].  See  Karlin  and  McGregor  (31). 
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Some  of  these  functions  are  spherical  functions  on  rank  one  spaces  over  reductive 
p-adic  groups.  See  Cartier  [18J  . 


cut  along  (-“,-1]  and 


lim  V>  <q  : q,x)  = ^(a;  x)  . 

q-i 

Another  interesting  case  is  q •+  -1.  Consider  the  case  8 = -|q|\  Set  q = -p. 


Then 


W (cos  9)  = 

-k|\q 


/ 2i0  , 
(e  ; q) 


, I | A 210  . 

( — I q I e :q> 


, 2i6  2.  , 2ie  2, 

(e  ; q l^fqe  ;q) 

, | jl  2i0  2.  , ill  2i0  2 

(- I q I e ;q  >„(-q|q|  « ;q  > 


, 2ie  2 , 2i0  2. 

(e  ; p ) (-pe  ;p  ) 


, 1+1  2i0  2 , 1 2i0  2. 

(p  e ;p  ) (— P e ; p ) 


I „ , -(1+1)  2 1+1  2i0.  „ , -(1-1)  2 X 2i0.| 

I i^o <p-  '•  p ' p e ’ ; P , -p  e )|  . 


As  p + 1 this  converges  to 


F 

1 0 


-<1  + l)/2 


2i9 


Set 


1F0 


(1-  D/2 


2i0 


2 ,,  2i0  (l+l)/2  2i0  (l-l)/2  2 

= (1  - e ) (1  + e ) 


(2-2  cos  20)  (X+1,/2(2  + 2 cos  20)  (X~1)/2 


,1+1,  . 2„, (1+11/2  1-1  2 (l-l)/2 

= 2 (sin  0)  2 (cos  0) 


,21,,  2,  (l+l)/2,  1 1-1 

= 2 (1  - x ) | x | 


(5.8) 


N* (cos  0)  = lim  C (cos  0;  -|q|*|q)  . 
n ^ , n 

q— -i 


Then 


(5.9) 

A change  of  variables  gives 

(5.10) 


/ N*(x)NA(x)  (1  - xV/2|x|X  ^dx  = 0,  m ^ n . 
^ m n 


(i,i-l, 

n'  (x)  = a P 2 2 (2x^  - 1)  , 

2m  m m 


(5.11) 


A A) 

v ' o 1 o ' 2 

n* tx)  " b xP  <2x  - 11  - 

2m+l  m m 


where  the  Jacobi  polynomial  P^°'^(x)  is  defined  by 

n 
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6.  The  continuous  q-Hermite  polynomials.  The  last  special  case  that  will  be  mentioned 


is  the  case  B = 0.  With  a different  normalization  they  will  be  called  the  continuous 


q-Hermite  polynomials 


The  recurrence  rel;  tion  becomes 


The  generating  function  (2.10)  becomes 


Rogers  (35)  and  later  Szego  [39)  and  Carlitz  [15) , [16)  studied  these  polynomials 


a slight  translation  to  correspond  to  those  of  Rogers  and  special  cases  of  formulas 


tions  are  polynomials  in  this  variable.  However  the  price  he  has  to  pay  of  losing 


the  nice  orthogonality  relation  is  so  high  that  we  will  use  the  notation  above 


The  q-extension  of  Mehler's  bilateral  sum  for  Hermite  polynomials  is 


i(9-*) 


See  [40,  Problem  23)  for  Mehler's  formula.  A beautiful  combinatorial  proof  of  Mehler's 


formula  has  been  given  by  Foata  [25)  . Formula  (6.5)  is  one  of  the  few  results  that 


have  been  found  for  q-Hermite  polynomials  that  has  not  been  extended  to  q-ultraspherical 


polynomials 


The  orthogonality  relation  (4.3)  becomes 


1 00  2tt  (q ; q) 

(6.6)  / H (x|q)H  (x|q)  TT  a-2(2x2-l)qk  + q2k)(l-xV1/2dx  = — — n-  6 

-1  " m k=0  “*»«>.  m'n 

so  (6.5)  is  the  analogue  of  Mehler’s  formula.  For  they  both  have  the  form 


l 


n=0 


Pn(x)P„<n)rn 
n n 

h 

n 


with 


b 

h = / [p  (x) ] 2w(x)dx  . 
n ' n 
a 

The  orthogonality  relation  (6.6)  was  given  by  Allaway  (2). 

The  reader  should  now  read  the  papers  of  Rogers  [34],  (35),  and  [36)  to  see  how 
he  used  the  q-Hermite  polynomials  to  obtain  the  Rogers-Ramanu jan  identities  and  many 
other  results. 

There  are  many  open  problems.  The  most  important  is  to  find  spaces  on  which 
these  polynomials  live,  presumably  as  spherical  functions.  Then  groups  acting  on 
these  spaces  need  to  be  used  to  derive  addition  formulas.  It  would  be  interesting  to 
find  a combinatorial  interpretation  of  the  q-Hermite  polynomials,  and  then  use  it  to 


give  a combinatorial  proof  of  (6.5).  Slepian  [38]  has  a very  important  multisum 
extension  of  Mehler’s  formula  and  Foata  and  Garsia  [26]  have  found  a combinatorial 


proof  of  this  formula.  A similar  formula  for  the  continuous  q-Hermite  polynomials 


should  be  obtained . 
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